#!/usr/bin/env python

'''
program plots Entropy
'''

import matplotlib.pyplot
import math
import numpy as np


T = [3.62, 3.618, 3.612, 3.61,  3.606, 3.6,   3.59,  3.58,  3.57,  3.56,  3.55,  3.53,  3.51,  3.49,  3.47,  3.45,  3.43,  3.41,  3.37,  3.34,  3.31,  3.275, 3.24,  3.20,  3.16,  3.11,  3.055, 3.0,   2.96,  2.91,  2.85,  2.79,  2.72,  2.65,  2.57,  2.48,  2.4,   2.32,  2.2,   2.1,   2,    1.9,   1.77,  1.61,  1.4,  1.2,    0.9,   0.4, 0]
#B = [0,    0.1,   0.2,   0.3,   0.4,   0.5,   0.6,   0.7,   0.8,   0.9,   1,     1.1,   1.2,   1.3,   1.4,   1.5,   1.6,   1.7,   1.8,   1.9,   2,     2.1,   2.2,   2.3,   2.4,   2.5,   2.6,   2.7,   2.8,   2.9,   3,     3.1,   3.2,   3.3,   3.4,   3.5,   3.6,   3.7,   3.8,   3.9,   4,    4.1,   4.2,   4.3,   4.4,   4.5,   4.6,  4.7, 4.76]
S = [0.31, 0.31,  0.308, 0.307, 0.307, 0.307, 0.304, 0.302, 0.301, 0.299, 0.298, 0.294, 0.289, 0.289, 0.283, 0.279, 0.276, 0.273, 0.265, 0.259, 0.255, 0.247, 0.241, 0.234, 0.227, 0.217, 0.207, 0.196, 0.192, 0.184, 0.174, 0.165, 0.155, 0.145, 0.136, 0.120, 0.111, 0.104, 0.089, 0.079, 0.072, 0.065, 0.051, 0.042, 0.026, 0.022, 0.008, 0, 0   ]

matplotlib.pyplot.xlabel("T")
matplotlib.pyplot.ylabel("S")
matplotlib.pyplot.title('Entropy along critical line \n $L = 27, $')
matplotlib.pyplot.plot(T,S)
matplotlib.pyplot.show()